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Abstract: When using the standard McCormick inequalities twice to convexity trilinear monomials, as is often 
the practice in modeling and software, there is a choice of which variables to group first. For the important 
case in which the domain is a nonnegative box, we calculate the volume of the resulting relaxation, as a 
function of the bounds defining the box. In this manner, we precisely quantify the strength of the different 
possible relaxations defined by all three groupings, in addition to the trilinear hull itself. As a by product, 
we characterize the best double-McCormick relaxation. 

We wish to emphasize that, in the context of spatial branch-and-bound for factorable formulations, our 
results do not only apply to variables in the input formulation. Our results apply to monomials that involve 
auxiliary variables as well. So, our results apply to the product of any three (possibly complicated) expressions 
in a formulation. 
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1. Introduction. Spatial branch-and-bound (sBB) (see [1, [2(1 [25]) for so-called factorable 
mathematical-optimization formulations (see |T5|) is the workhorse general-purpose algorithm in 
the area of global optimization. It works by using additional variables to reformulate every function 
of the formulation as a (labeled) directed acyclic graph (DAG). Root nodes can be very complicated 
functions, and leaves are variables that appear in the input formulation, each labeled with its 
interval domain. Intermediate nodes are labeled with auxiliary variables together with operators 
from a small dictionary of basic functions of few (often one or two) variables. Also, we have a 
method for convexifying the graph of each dictionary function. sBB algorithms work by composing 
convex relaxations of the dictionary functions, according to the DAG, to get relaxations of the root 
functions. Bounds on the leaves propagate to other nodes and conversely. Branching (subdividing 
the domain interval of a variable) creates subproblems, which are treated recursively. Objective 
bounds for subproblems are appropriately combined to achieve a global-optimization algorithm. 

Much of the research on sBB has focused on developing tight convexifications for basic functions 
of few variables (many references can be found in [5j). Other research has focused on how bounds 
can be efficiently propagated and how branching can be judiciously be carried out (see J3] , for 
example). From the viewpoint of good convexifications, much less attention has been paid to how 
the DAGs are created, but this can have a strong impact on the quality of the resulting convex 
relaxation of the input formulation; see P2H31 (22 [29] for some key papers with other viewpoints 
concerning constructing DAGs. 

For basic multilinear monomials f(x 1 ,...,x n ) := X\ - ■ ■ x n , with Xi £ [a*,^], there is already a 
lot of flexibility which can have a significant impact on the overall convexification of the graph of 
f(x i,..., x n ) := X\ - ■ ■ x n on the box domain [aq, &i] x • • • x [ a n , b n \. For n = 2, we have the classic 
McCormick inequalities (see m ), which is simply the tetrahedron that is the convex hull of 
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the points (/, x±,x 2 ) '■= (a 1 a 2 , cq, a 2 ), (a 1 b 2 ,a 1 ,b 2 ), (b 1 a 2 ,bi,a 2 ), (bib 2 ,bi,b 2 ). The inequalities can be 
derived from the four inequalities 


{xi - a 1 )(x 2 - a 2 ) > 0, (xi - ai)(b 2 - x 2 ) > 0, 

(&i - x 1 )(x 2 - a 2 ) > 0, (6i - Xi){b 2 - x 2 ) > 0, 

by multiplying out and then replacing all occurrences of X\X 2 by the variable /. 

For general n, there are 2” points to consider (i.e., all choices of each variable at a bound), and 
the inequality descriptions in the space of (/, xi,... , x„) G M" +1 get rather complicated (even for 
n = 3; see pa ina). It is frequent practice, both in modeling and software, to repeatedly use the 
McCormick inequalities when n > 2. Already the trilinear case, n = 3, is an interesting one for 
analysis. Here, we have three choices, which can be thought of as / = (xix 2 )x 3 , / = (xix 3 )x 2 and 
/ = (x 2 x 3 )xi. Because the domain of each variable is its own interval \a t , h ( ], the grouping can affect 
the quality of the convexification. In what follows, we analytically quantify the quality of these 
different convexification possibilities, in addition to the trilinear hull itself. 

Our results are not just relevant to trilinear monomials in formulations. With the sBB approach 
for factorable formulations, our results are relevant whenever three quantities are multiplied. That 
is, as an expression DAG is created and auxiliary variables are introduced, a trilinear monomial will 
arise whenever three quantities (which can be complicated functions themselves) are multiplied. 

In what follows, we use (n + l)-dimensional volume to compare different natural convexifications 
of graphs of functions of n variables on the box domain [ai, bf\ x • ■ • x [o„, b n \. We present a complete 
analytic analysis of the case of n = 3, for all choices of 0 < a, < bi. It is perhaps surprising that this 
can be carried out, and probably less surprising that the analysis is quite complicated. 

Volume as a measure for comparing relaxations was first proposed in HQ]; also see [8] and |28j. 
In fact, the practical use of volume as a measure for comparing relaxations in the context of 
nonlinear mixed-integer optimization, foreshadowed by B33. was later validated computationally 
for a nonlinear version of the uncapactitated facility-location problem (see 0). Specifically, using 
volume calculations, a main mathematical result of m is that weak formulations of facility-location 
problems are very close to strong formulations when the number of facilities is small compared to 
the number of customers. Then [9] showed that in this scenario, with a convex objective function, 
the weak formulation computationally out performs the strong formulation in the context of branch- 
and-bound. The emphasis in mmm was not on sBB nor on low-dimensional functions. Because 
those results pertained to varying dimension and related asymptotics, exactly how volumes are 
compared and scaled was important (in particular, see [10] which defines the “idealized radial 
distance”). Because we now focus on low-dimensional polytopes, the exact manner of comparison 
and scaling is much less relevant. Using volume as a measure corresponds to a uniform distribution 
of the optimal solution across a relaxation. This is justified in the context of nonlinear optimization 
if we want a measure that is robust across all formulations. One can well find situations where 
the volume measure is misleading. It would not make sense for evaluating polyhedral relaxations 
of the integer points in a polytope, if we were only concerned with linear objectives — in such 
a case, solutions are concentrated on the boundary and there are better measures available (see 
urn- But if we are interested in a mathematically-tractable measure that robustly makes sense in 
the context of global optimization, volume is quite natural. 

Motivated by two well-studied applications (the Molecular Distance Geometry Problem and the 
Hartree-Fock Problem ), pfj first proposed volume in the context of sBB and monomials, but they 
leapfrogged to the case of n = 4 and took a mostly experimental approach. They demonstrated 
that there can be a significant difference in performance depending on grouping, and they offered 
some guidance based on computational experiments. At the time of that work, it appeared that 
developing precise formulae for volumes relevant to repeated McCormick was not tractable. With 
our present work on n = 3, it now seems possible that the case of n = 4 could be carried out (see fj9] 
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for an idea concerning how our results for n = 3 could already be applied in practice to the n = 4 
case). 

There has been considerable research on multilinear monomials and generalizations in the context 
of global optimization, notably pn 0 m, ed 0 M- Our work adds to that literature. 

In |2j we define the polytopes that we work with. In t|3j we discuss the various alternatives for 
working with triple products. In Q we present our main results and their consequences. In 
we present our proofs. In fj9j we describe future directions for investigation. <10 
contains technical lemmas and calculations. 


an appendix, 


2. Double McCormick. When using the double-McCormick technique to convexity trilinear 
monomials, a modeling/algorithmic choice is involved: we must choose to which pair of variables we 
apply the first iteration of McCormick. For the variables Xi 6 [a,, bi], i = 1,2,3, let Oi := afbjbk) + 
bfdjdk)- Then we can label the variables such that 0\ < 0 2 < O 3 . In this manner, we can assume 
that 

dib 2 b 3 + &ia 2 a 3 < 6ia 2 6 3 + aib 2 a 3 < bib 2 a 3 + a 3 d 2 b 3 . (II) 

Given the trilinear monomial / := x 3 x 2 x 3 , there are three choices of convexifications depending 
on the bilinear sub-monomial we convexity first. We could first group Xi and x 2 and convexity 
w = X\X 2 ] after this, we are left with the monomial / = wx 3 which we can also convexity using 
McCormick. Alternatively, we could first group variables X\ and x 3 , or variables x 2 and x 3 . 

2.1. Convexification. To see how to perform these convexifications in general, we show the 
double-McCormick convexification that first groups the variables x t and Xj. Therefore we have 
/ = XiXjX k and we let w t j = XiXj so / = WijX k . 

Convexifying w t j = XiX 3 we obtain the inequalities: 

— djXi — diXj + didj> 0, 

—+ bjXi + d^j — dibj > 0 , 

—w^ + djXi + biXj — bidj > 0 , 

w^ — bjXi — b^j + bibj > 0. 

Convexifying / = w^x k we obtain the inequalities: 

f dfcWij didjX k -f“ didjd k T 0; 

f T b k Wij T didjX k didjb k T 0; 

-/ + d k Wij + b.ibjX k - bibjd k > 0 , 

/ - b k Wij - bibjX k + bibjb k > 0 . 


Using Fourier-Motzkin elimination, we then eliminate the variable to obtain the following 
system in our original variables /, Xi , Xj and x k . 


Xi 

di 



> 

0, 

(1) 

Xj 

—dj 



> 

0, 

(2) 

f 

djdkXi 

didkXj 

- didjXk + 2 CL^CLjCL^ 

> 

0, 

(3) 

f 

djbkXi 

dib^x j 

bibjX k + didjbk + bibjbk 

> 

0, 

( 4 ) 

— 

Xj + bj 



> 

0, 

( 5 ) 

— 

Xi + bi 



> 

0, 

(6) 

f 

bjd k Xi 

bidkXj 

didjXk “l - didjdk + bibjdk 

> 

0, 

(7) 

f 

~bjb k Xi 

- bib k Xj - 

bibjXk + 2bibjb k 

> 

0, 

( 8 ) 

— 

f + bjbkXi + dib k x j 

; + didjXk didjbk dibjbk 

> 

0, 

(9) 

- 

/ + djb k Xi + bihxj 

; + didjXk didjbk bidjbk 

> 

0, 

(10) 
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-x k + b k 

>0, 

(ID 

f + bjd k Xi “h dia k Xj + bibjX k OjibjCt k bibjd k 

>0, 

(12) 

f + djd k Xi + bid k Xj + bibjX k bi&jCL k bibjd k 

>0, 

(13) 

k ctk 

>0, 

(14) 

Q>i(LjX k 

>o, 

(15) 

- / + bibjX k 

>0. 

(16) 


16 


'OfcQ + did). , and 


It is easy to see that the inequalities 15 and 16 are redundant: 15 is a 

+ ( fl2| ). 

We use the following notation in what follows. For i= 1, 2,3, system i is defined to be the system 
of inequalities obtained by first grouping the pair of variables Xj and Xk, with j and k different 
from i. Vi is defined to be the solution set of this system. 


2.2. Hull. As we noted earlier, a convex-hull representation for trilinear monomials is known. 
From m, for any labeling that satisfies [TI] (or even just the first inequality of[TI]), this inequality 
system which we refer to as H is: 


/ - d 2 a 3 x i - a 1 d 3 x 2 - d x d 2 x 3 + 2d x d 2 a 3 > 0, (17) 

/ - b 2 b 3 x-i - bib 3 x 2 - b l b 2 x 3 + 2bib 2 b 3 > 0, (18) 

/ — d 2 b 3 X\ — d\b 3 x 2 — bid 2 x 3 + d\d 2 b 3 + b\d 2 b 3 > 0, (1-9) 

/ — b 2 d 3 X\ — b\d 3 x 2 — d\b 2 x 3 + b 3 b 2 d 3 + dib 2 a 3 > 0, (20) 

/ — bf-a 1 ~ b\d 3 x 2 — bid 2 x 3 + ^ bl 1 -^ + b\b 2 a 3 + b\d 2 b 3 — ai6 2 6 3 ^ > 0, (21) 

/ ~~ a ] n -b 1 x i ~ a ib 3 x 2 — dib 2 x 3 + ^ df-bx + °i a 2&3 + d\b 2 d 3 — bid 2 d 3 j > 0, (22) 

—f d 2 d 3 Xi + b 3 d 3 x 2 + b\b 2 x 3 — b 3 b 2 d 3 — b±d 2 d 3 > 0, (23) 

—/ + b 2 d 3 X\ + d 1 d 3 x 2 + b\b 2 x 3 — bib 2 d 3 — dib 2 a 3 > 0, (24) 

—/ + d 2 a 3 X\ + bib 3 x 2 + b\d 2 x 3 — b 3 d 2 b 3 — b 3 d 2 d 3 > 0, (25) 

—/ + b 2 b 3 Xi + aia 3 x 2 + dib 2 x 3 — d\b 2 b 3 — d\b 2 a 3 > 0, (26) 

—/ + d 2 b 3 X\ + bib 3 x 2 + d 3 d 2 x 3 — bid 2 b 3 — a 3 d 2 b 3 > 0, (27) 

—/ + b 2 b 3 x i + d\b 3 x 2 + d 3 d 2 x 3 — d\b 2 b 3 — a\d 2 b 3 > 0, (28) 

X\ — d\ > 0, (29) 

— X\ + bi > 0, (30) 

x 2 -a 2 > 0, (31) 

—x 2 + b 2 > 0, (82) 

x 3 -a 3 > 0, (33) 

—x 3 + b 3 > 0, (34) 


where = bib 2 d 3 - a 3 b 2 b 3 - bid 2 d 3 + bid 2 b 3 and ?? 2 = cq d 2 b 3 - b x d 2 d 3 - ai6 2 ^ 3 + ai^a 3 . 


We refer to the polytope defined as the feasible set of system H as Vh- The extreme points of 
Vh are the 8 points that correspond to the 2 3 = 8 choices of each x-variable at its upper or lower 
bound. We label these 8 points (all of the form [f = xix 2 x 3 ,x 1 ,x 2 ,x 3 \) as follows: 



Oia 2 a 3 


a 1 a 2 b 3 


ai b 2 a 3 

2 

ai 

3 

ai 

4 

ai 

, v := 

a 2 

, v := 

a 2 

, v := 

b 2 


°3 


b 3 


a 3 


d\b 2 b 3 


b\b 2 b 3 


b\b 2 a 3 


b\a 2 b 3 

Oi 

6 . 

bi 

7 . 

hi 

8 . 

h 

b 2 

, v := 

^2 

, v := 

b 2 

, v := 

a 2 

b 3 


&3 


a 3 


b 3 
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Each alternative polyhedral convexification leads to a different system of inequalities (system i, 
i = 1,2,3) and therefore a different polytope (T l , i = 1,2,3) in M 4 — all three contain the convex 
hull of the solution set of our original trilinear monomial (on the box domain), i.e. Vr. 

To establish if one of these three convexifications is better than another, we need to be able 
to compare these polytopes in a quantifiable manner. We take the (4-dimensional) volume as our 
measure, with the idea that a smaller volume corresponds to a tighter convexification. 

For trilinear monomials with domain being a box (in the nonnegative orthant), we derive exact 
expressions for the (4-dimensional) volume for the convex hull of the set of solutions and also for 
each of the three possible double-McCormick convexifications. These volumes are in terms of six 
parameters (the upper and lower bounds on each of the three variables) and are rather complicated. 
By comparing the volume expressions, we are able to draw conclusions regarding the optimal way 
to perform double McCormick for trilinear monomials. 


3. Alternatives. In practice, there are many possibilities for handling each product of three 
terms encountered in a formulation. A good choice, which may well be different for different triple 
products in the same formulation, ultimately depends on trading off the tightness of a relaxation 
with the overhead in working with it. For clarity, in the remainder of this section, we focus on 
different possible treatments of / = x\x 2 x 3 . 

One possibility is to use the full trilinear hull Vr. This representation has the benefit of using no 
auxiliary variables. Another possibility to use the convex-hull representation (see [5], for example), 
writing / = )T® = i A/C, with )T* =I A^ = 1, A j > 0, for j = 1,2,..., 8. This formulation has the draw¬ 
back of utilizing eight auxiliary variables. But noticing that there are 5 linear equations, we can 
really reduce to three auxiliary variables. In fact, there is a very structured way to do this, where 
none of the A j variables are employed at all, and rather we introduce three auxiliary variables w 12 , 
Wi 3 and w 2 3, which represent the products XiX 2 , x 3 x 3 and x 2 x 3 , respectively. A strong advantage 
of this last approach is when terms x 3 x 2 , x 3 x 3 and x 2 x 3 are also in the model under consideration. 
We wish to emphasize that projecting any of these convex-hull representations (reduced or not) 
down to the space of (f,x 1 ,x 2 ,x 3 ) yields again Vr, and so all of these representations have the 
same bounding power. 

We are advocating the consideration of double-McCormick relaxations as an alternative when 
warranted. We have identified the best among the double McCormicks and quantified the error in 
using it in preference to Vr (and, ipso facto, with any convex-hull or reduced convex hull represen¬ 
tation). A double-McCormick relaxation involves only one auxiliary variable (and 8 inequalities). 
This can be particularly attractive when this particular auxiliary variable already appears in the 
model under consideration. Alternatively, especially when this particular auxiliary variable does 
not appear in the formulation, we can use the formulation with zero auxiliary variables QH- 
Recently (see E3), we have computationally validated such an approach in the context of “box 
cubic programs” 


mm 

i6i" 


,n 


^ ^ Clijk J'i J'j %k • ?Ci £ fail bi\ , i 1 , 2 ,., 

_ {iJM 

In this type of problem, we can apply (T][l4) independently for each trinomial, with no auxiliary 
variables at all, choosing the best double-McCormick for each trinomial, whenever the associated 
volume is close to the volume for Vr. We have documented that this can happen quite a lot, and so it 
is a viable approach. It is important to emphasize that some of the negative experience with double 
McCormick is related to choosing the wrong one. Indeed, our mathematical and computational 
results indicate that there are many situations where: (i) the worst double McCormick is quite bad 
compared to the best one, and (ii) the best one is only slightly worse than Vr (and its convex-hull 
representations). 
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Besides any prescriptive use of double-McCormick relaxations, our results can simply be seen 
as quantifying the bounding advantage given by Vh and the various convex-hull representations 
(reduced or not) as compared to each of the possible double McCormick relaxations. 

In some global-optimization software (e.g., BARON and ANTIGONE) the complicated inequality 
description of the trilinear hull is explicitly used. In other global-optimization software (e.g., 
COUENNE and SCIP) and as a technique at the formulation level, repeated McCormick is used for 
the trilinear case. It is by no means clear that either approach should be followed all of the time 
(though this currently seems to be the case), because of the solution-time tradeoff in using more 
complicated but stronger convexihcations. This effect can be especially pronounced in the case of 
nonlinear optimization where solutions may not be on the boundary (see P, for example). By 
quantifying the quality of different convexihcations, we offer (i) firm and actionable means for 
deciding between them at run time and, (ii) some explanation for differing behavior of sBB software 
under different scenarios. 

Finally, we note that the double-McCormick approach is often applied at the modeling level 
(see m , for example). In particular, our results are highly relevant to modelers who simply use 
global-optimization software, often through a modeling language. An uniformed modeler can defeat 
clever software. In such a case, it is very useful for the user to know which double McCormick to 
use, because a bad one may negatively affect sBB performance, and all sBB software that we know 
will not capture implicit triple products in formulations. 

4. Theorems. First we define the following twelve points in M 4 , where j :=i+ 1 (mod 3) and 
k := i + 2 (mod 3): 


where: 


01 

0- 

oj 


rail 


re?l 


refl 


r^i 


rflJi 


rail 

el 

10 

e\ 

11 

Of 

12 

el 

9 

bi 

10 

o\ 

0,1 

,V1 

bi 

,vi ■— 

b2 

,Vj ■— 

02 

,v 2 

el 

,v 2 

el 

^3 


a 3 


03 


b 3 


03 


^3 



\ell 


\el] 


COCO 


p!l 


co-a 



rflsi 

11 

Oi 

12 

bi 

9 

b 1 

10 

Oi 

11 

ffli 

12 


bi 

v 2 := 

el 

bi 

,v 2 := 

el 

a 3 

,V 3 ■■= 

a 2 

et 

,V 3 ■■= 

b2 

el 

,^3 := 

b2 

el 

,V 3 


a 2 

el 


djdjd k + 


djdjd k T 


dj(b k a k ){bibjb k djdjd k ) 

bjb k Qj CL k 

a k (bj dj')(J)ibjb k diCLjd k ) 

bj b k dj d k 

bjd k (cijbjb k didjb k bidjd k T b^djb k ) 

bj b k dj d k 

djb k {bibjd k bidjdk (i) bj o k O:, bj b k ) 
bj b k dj d k 


Qi = dj + 

Of = a i + 

Oi = di + 

&i = dj + 


dj(bi df)(b k d k ^j 

bj b k dj d k 
d k (bj-dj)(bi-di) 

bj b k dj d k 
bjjbj- dj)(b k - d k ) 

bj b k dj d k 
b k {b j ~dj){b l -d i ) 

bj b k dj d k 


Next, we state our main results. 


Theorem 4.1. 

Vol- p H = (bj — d\)(b 2 — 02 X ^3 ~ R 3 ) x 

( 61 ( 562^3 — ®2^3 — & 2 O 3 — 30203 ) + 01(50203 — 6203 — 02&3 — 36263 )) /24. 

Theorem 4.2. The set of extreme points of V\ is {v 1 ,... ,u 8 } U {u®,... ,uj 2 }. Moreover, 

Vol-p 1 = Vol -p H +( 6 i — Oi )(62 — 02 X (^3 — dff x 

3(636203 — 036203 + 630263 — 030263 ) + 2(036363 — 630203 ) 

24(6263 — 0203 ) 
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Theorem 4.3. The set of extreme points ofV 2 is {u 1 , ... ,u 8 } U {u®, ..., vf 2 }. Moreover, 


Vol-p 2 — Vol -p H + 


(61 - Qi)(6 2 - a 2 ) 2 (^3 - Q 3) 2 (5(aib 1 b 3 - ai 6 ia 3 ) + 3(6 2 a 3 - a 2 6 3 )) 

24(6163 - aia 3 ) 


Theorem 4.4. The set of extreme points ofV 3 is {u 1 ,... ,u 8 } U {uf,..., u] 2 }. Moreover, 


Volp 3 — Vol -p H + 


(61 - ai)( 6 2 - a 2 ) 2 ( 6 3 - a 3 ) 2 ( 5(016362 - ai 6 ia 2 ) + 3(6® a 2 - a 2 6 2 )) 


24(6i6 2 — Oid 2 ) 

Our proofs in §^5]^8] all assume that ai,a 2 ,a 3 > 0. Next, we briefly explain why the theorems 
hold even when any of the a, are zero. Taking the convex hull of a compact set is continuous 
(even 1-Lipschitz) in the Hausdorff metric (see [ 251 p. 51]). The volume functional is continuous 
(with respect to the Hausdorff metric) on the set K n of convex bodies in R” (see {23], Theorem 
1.8.20; p. 68 ]). If two sets of m points in R" are close as vectors in R mn , then they are also close 
in the Hausdorff metric. Therefore, the volume of the convex hull of a set of m points in R ra is a 
continuous function of the coordinates of the points. Also, the coordinates of the extreme points 
of our polytopes are all continuous functions (of the six parameters) at a, = 0. Finally, we note 
that the volume formulae that we derive are continuous functions (of the six parameters) at a* = 0 . 
Therefore, those formulae are also correct when some a* = 0. We do note that we can also modify 
our constructions to handle these cases where some of the a* are zero, but our continuity argument 
is much shorter. 


Corollary 4.1. For all values of the parameters ai, 61 ,a 2 , 6 2 ,a 3 , 6 3 , meeting the conditions 
(foj) . we have: Vo\-p H < Vofp 3 < Volp 2 < Vol- Pl . 

From this we can see that with the variables ordered according to their upper and lower bounds 
per ([H]), the smallest volume will always be obtained by using system 3 (i.e., first grouping variables 
Xi and x 2 ). In addition, for different values of the upper and lower bounds, we can precisely quantify 
the difference in volume of the alternative convexifications. 

Moreover, by substituting oq = a 2 = 0 and 61 = b 2 = 1 into the conditions ([H]), we can easily see 
the following corollary relevant to mixed-integer nonlinear optimization. 

Corollary 4.2. In the special case of two binary variables and one continuous variable, first 
grouping the two binary variables gives the convexification with the smallest volume. 

In this special case, we only have two parameters a 3 and 6 3 and the volume formulae simplify 
considerably. In particular, for this special case, V 3 is equivalent to Vh, and V\ and V 2 are equiva¬ 
lent. We compute the difference in volume between the two distinct choices of convexification and, 
in Figure [lj plot this expression as the parameters vary (satisfying 0 < a 3 < 6 3 ). The following is 
easy to establish. 

Corollary 4.3. As a 3 and b 3 increase, the difference in volumes ofV 3 and V 2 (orV 3 ) becomes 
arbitrarily large. Additionally, for a fixed b 3 , the greatest difference in volume occurs when a 3 = b 3 /3. 

Finally we note that in the special case in which a 3 = a 2 = a 3 = 0, each convexification reduces to 
the convex hull, which is a result of [2Tlj . So in this case, any double-McCormick convexification has 
the power of the more-complicated inequality description of the convex hull. In fact, viewed this 
way, our results provide a quantified generalization of this result of m- We do wish to emphasize 
that because our results do not just apply to trilinear monomials on the formulation variables, but 
may well involve auxiliary variables, the case of non-zero lower bounds is very relevant. 





Speakman and Lee: Quantifying Double McCormick 


Figure 1. Graph of difference in volume 


( 3a 3 (b 3 -a 3 ) 2 

V 24t>3 


vs. parameter values 



5. Proof of Thm. |4.1[ We compute the volume of Vh by constructing a triangulation. See 
Figure [ 2 ] for a diagram of the 8 extreme points of Vh- Note that v 2 , which has all of the variables 
at their lower bounds, is at the bottom of the “inner cube”, and v 6 , which has all of the variables 
at their upper bounds, is at the top of the “outer cube”. 

We use the fact that the volume of an n-simplex in K” with vertices (z°,..., z n ) is: 

| det(^ 1 — z° z 2 -z° ... z n -z°)\/n\. 

The volume of the 4-simplex with extreme points v 1 ,v 2 ,v 4 ,v 5 and v 6 , which we define as S ■= 
convju 1 , v 2 , u 4 , v 5 , u 6 }, is 


(61 - ai) 2 ( 6 2 - a 2 )(b 3 - a 3 )(b 2 b 3 - a 2 a 3 )/24. 

A 4-simplex has 5 facets, each of which is a 3-simplex and is described by the hyperplane through 
a choice of 4 extreme points. To determine the facet-describing inequalities, we compute each 
hyperplane and then check the final point to obtain the direction of the inequality. The 5 facets of 
S are described as follows: 

F 1 (hyperplane through points v 1 , v 2 , u 4 , v 6 ): 

, , , (aia 2 a 3 - ai b 2 a 3 - bia 2 a 3 + M 2 6 3 ) 

—J + o, 2 a 3 X\ + a 3 a 3 x 2 -| - — -, - x 3 

( b 3 — a 3 ) 

(a\a 2 a 3 b 3 - a x b 2 a 2 3 - b^aj + 6 i 6 2 a 3 6 3 ) ^ 

(b 3 - o 3 ) 

F 2 (hyperplane through points v 1 ,v 2 ,v A ,v 5 ): 

f — a 2 a 3 Xi — a 3 a 3 x 2 — ®i b 2 x 3 + CLia, 2 a 3 + cii 6 2 o 3 ^ 0 
F 3 (hyperplane through points v 1 , v 2 , v 5 , v 6 ): 


( b 3 - a 3 )x 2 - (b 2 - a 2 )x 3 + b 2 a 3 - a 2 b 3 > 0 
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Figure 2. Visual representation of simplex, 5, with extreme points v 1 ,v 2 ,v 4 ,v 5 and v 6 , plus the other 3 convex hull 
extreme points v 3 , v 7 and v 8 . 



F 4 (hyperplane through points u 1 ,'i; 4 , , u 5 ,u 6 ): 


/ - 6 2 6 3 Xi - 


— 0162X3 + 


( aib 2 a 3 - aib 2 b 3 - feia 2 Q 3 + M263) 

(b 2 -a 2 ) 

(— a,\CL 2 b 2 b 3 + (iib^cis — bia 2 b 2 a 3 + 616263) 


(62 - a 2 ) 


> 0 


-F 5 (hyperplane through points v 2 , v 4 , v 5 , v 6 ): 

—/ + 6263X1 + 0-103X2 + 0162X3 — O162O3 — 036263 > 0 


If a hyperplane if intersects a polytope P on a facet F, then H + (resp., H ~) denotes the half¬ 
space determined by H that contains (does not contain) P. If a point w is not in H but in H + 
(resp., H~), then w is beneath ( beyond ) F (see [6j p. 78]). 

We now compute the volume of conv(5U {u 8 }). To obtain the additional volume of this polytope 
compared with S we sum the volume of conv({u 8 } U F) for each facet, F, of S such that v 8 is 
beyond that facet. To do this we first check each of the 5 facets to determine if v 8 is beneath or 
beyond that facet. To do this we substitute v 8 into the relevant inequality, and if the result is 
negative then v 8 lies beyond that facet. 

It is easy to check that v 8 satisfies F 1 and F 5 and violates F 3 . Using Lemma 10.1 , we also check 
that v 8 satisfies both F 2 and F 4 . From this we have that v 8 is beyond one facet, F 3 . Therefore we 
need to calculate the volume of conv(F 3 U {u 8 }) = conv{u 1 ,r; 2 ,u 5 ,?; 6 ,u 8 }, this a a 4-simplex with 
volume: 


(61 - ai) 2 (6 2 - a 2 )(6 3 - a 3 )(6 2 6 3 - a 2 a 3 )/24. 


We now have a new polytope which is convju 1 , v 2 , v 4 , v 5 , v 6 , u 8 } = conv(5 U {u 8 }). We refer to 
this polytope as Q. The volume of Q is given by the sum of the volumes of the two simplices 
we have computed thus far. The facets of Q are the facets of the original simplex without F 3 , 
along with the facets of the 4-simplex: conv(F 3 U {u 8 }) (again not including F 3 itself). A facet of 
conv(F 3 U {u 8 }) is supported by a hyperplane through a choice of 4 of the 5 extreme points (points 
v 1 ,v 2 ,v 5 ,v 8 and v 8 ). As before, to determine these facet inequalities we compute each hyperplane 
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and then check the final point to obtain the direction of the inequality (note that we exclude the 
choice v l ,v 2 ,v 5 ,v 6 because this corresponds to F 3 ). The 4 facets are described below: 

F 6 (plane through points v 1 , v 2 , v 5 , v 8 ): 


(—Oia 2 a 3 + aib 2 b 3 + &ia 2 a 3 - M 2 & 3 ) 

J - a 2 a 3 x i--r- x 2 

(o 2 — a 2 ) 

(—a-ia 2 a 3 + cli a 2 b 2 b 3 — b\a\b 3 + b\a 2 b 2 a 3 ) 

— t>ia 2 x 3 H- — -,- 

[b 2 - a 2 ) 

F‘ (plane through points v 1 ,v 5 ,v 6 ,v 8 ): 

f - b 2 b 3 x i - b 3 b 3 x 2 - b 1 a 2 x 3 + b ± a 2 b 3 + &iM 3 > 0 


> 0 


F 8 (plane through points v 2 ,v 5 ,v 6 ,v 8 ): 

-f + b 2 b 3 x i + bib 3 x 2 + 


(— a 3 a 2 a 3 + a 3 b 2 b 3 + b\a 2 b 3 — bib 2 b 3 ) 


(b 3 — a 3 ) 

(—aia 2 a 3 b 3 + a 3 b- 2 b 2 + b\a 2 b\ — bib 2 a 3 b 3 ) ^ 

(b^f) - 

F 9 (plane through points v 1 ,v 2 ,v 6 ,v 8 ): 

—/ + a 2 a 3 x i + bib 3 x 2 + bia 2 x 3 — b\a 2 a 3 — bia 2 b 3 > 0 


~x 3 


The facets of Q = convjn 1 , v 2 ,v 4 , v 5 , v 6 , n 8 } are therefore F 1 , F 2 , F 4 , F 5 , F 6 , F 7 , F 8 and F 9 . 

To obtain the entire volume of Vh w e need to consider two further extreme points: v 3 and v 7 . 
It would be convenient to add these points separately; i.e., compute the additional volume each 
produces when added to Q, and sum the results. As the following lemma shows, this will give the 
correct volume if the intersection of the line segment between these points and Q is not empty. 

Lemma 5.1. Let P be a convex polytope and let w\ and w 2 be points not in P. Let L(wi, vj 2 ) be 
the line segment between w 3 andw 2 . If L(yj\, w 2 )nP / 0, then conv(P, w 1 )L)conv(P, w 2 ) is convex. 
Moreover, in this case, conv(P,Wi,w 2 ) = conv(P,w 1 )L) conv(P,w 2 ). 

Proof. First we show that conv{P,Wi) U conv(P,w 2 ) is convex. If we show that L(w 3 ,w 2 ) is 
completely contained in conv(P,Wi)L)conv(P,w 2 ), then we will be done. Choose z £ L{wi,w 2 )C]P. 
Now consider L{w\,z). Because z € P, this whole line segment must be in conv{P,W\). Similarly 
consider L(z,w 2 ); this whole line segment must be contained in conv(P, w 2 ). Therefore the whole 
line segment L(w 3 ,w 2 ) must be contained in conv(P,Wi) U conv(P,w 2 ) and therefore this set is 
convex. 

Next, we demonstrate that conv(P,Wi,w 2 ) = conv(P,w 1 ) U conv(P,w 2 ). First, choose y £ 
conv{P,Wi ) U conv(P,w 2 ); therefore y £ conv{P,Wi) or y £ conv{P,w 2 ) (or both); in either case 
it is clear that y £ conv(P,Wi,w 2 ). In the other direction, choose y £ conv(P, wi,w 2 )] therefore 
y can be written as a convex combination of the extreme points of P and w\ and w 2 . Because 
conv(P,Wi ) U conv(P,w 2 ) is convex, this means y £ conv{P,wf) U conv{P,w 2 ). Therefore the sets 
are equal as required. □ 

We refer to the midpoint of the line between w 3 and w 2 as M{wi,w 2 ). To show that the inter¬ 
section of L(v 3 ,v 7 ) and Q is non-empty, consider the midpoint 

M(y 3 V 7 ) = [ a l a 2 b 3+bib2a3 ai+b 1 a2+&2 b 3+ a 3 j 

We show that this point satisfies each of the inequalities of Q by substituting into each inequality 
and checking the result. By showing that each resulting quantity is nonnegative, we conclude 
that the midpoint intersects Q. It is easy to see that the midpoint M(v 3 ,v 7 ) satisfies F 1 , F 5 , F 8 
and F 9 . Using Lemma |10.1[ we also check that M{v 3 ,v 7 ) satisfies F 2 , F 4 , F 6 and F 7 . Therefore 
conv(QU {n 3 }) U conv(Q U {T 7 }) = conv(QU {v 3 } U {c 7 }) = Vh- 
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5.1. Computing the (additional) volume of conv(Q U {n^}). We now compute the 
additional volume of conv(Qu{u 3 }) compared to the volume of Q. To obtain this, we sum the 
volumes of conv({u 3 } U F) for each facet, F, of Q such that v 3 is beyond that facet. We substitute 
v 3 into each relevant inequality, and if the result is negative then v 3 lies beyond that facet. It is easy 
to see that v 3 satisfies F 5 , F 7 and F 9 and violates F 2 , F 6 and F 8 . It can then be checked that v 3 


satisfies F 1 using Lemma 10.4 (with A = b 2 , B = a 2 , C = ( 6 i 6 3 — aia 3 ), D = (2aia 3 — ai 6 3 — feia 3 )). We 


also check that v 3 satisfies F 4 using Lemma 10.4 (with A = b 2 , B = a 2 , C = (a 3 a 3 — 2ai6 3 + 6 1 6 3 ),.D = 


(aib 3 — &ia 3 )) and Lemma 10.1 


From this we know that v 3 is beyond F 2 , F 6 and F 18 ; therefore we need to compute the volume 
of the convex hulls of v 3 with each of these facets. 

The polytope conv(T 2 U {w 3 }) = convju 1 , v 2 , v 4 , v 5 , u 3 } is a 4-simplex with volume: 

ai( 6 i - ai )(&2 - a 2 ) 2 (b 3 - a 3 ) 2 /24. 

The polytope conv(F 6 U {u 3 }) = convju 1 , v 2 , v 5 , v 8 , u 3 } is a 4-simplex with volume: 

02(61 — Oi ) 2 (62 — a 2 )( 6 3 — a 3 ) 2 /24. 

The polytope conv(F 8 U {u 3 }) = conv{v 2 ,v 5 ,v 6 ,v s ,v 3 } is a 4-simplex with volume: 

6 3 ( 6 i — Oi) 2 (6 2 — a 2 ) 2 (b 3 — o 3 )/24. 

5.2. Computing the (additional) volume of conv(Q U {u^}). We now compute the 
additional volume of conv(Q U {u 7 }) compared to the volume of Q. To obtain this, we sum the 
volumes of conv({u 7 } U F) for each facet, F, of Q such that v 7 is beyond that facet. We substitute 
v 7 into each relevant inequality, and if the result is negative then v 7 lies beyond that facet. It is 
easy to see that v 7 satisfies F 2 , F 5 and F 9 and violates F 1 , F A and F 7 . It can then be checked that 


v 7 satisfies F 6 using Lemma 10.4 (with A = b 2 ,B = a 2 ,C = ( 6 3 a 3 — a 3 b 3 ),D = (a 3 a 3 — 2b 1 a 3 + bib 3 )) 
We also check that v 7 satisfies F 8 using Lemma [lO .4 (with A = 6 2 , B = a 2 ,C = 


10.1 


and Lemma 

(26i 6 3 - ai 6 3 - 61 a 3 ), D = (a ± a 3 - &i6 3 )). 

From this we know that v 7 is beyond F 1 , F 4 and F 7 , therefore we need to compute the volume 
of the convex hulls of v 7 with each of these facets. 

The polytope conv(T 1 U {u 7 }) = convju 1 ,?; 2 ,!; 4 ,?; 6 ,!: 7 } is a 4-simplex with volume: 

0 - 3(61 — Oi ) 2 (62 — a 2 ) 2 ( 6 3 — a 3 )/24. 

The polytope conv(T 4 U {w 7 }) = convlu 1 ,?; 4 ,!; 5 ,!; 6 ,!; 7 } is a 4-simplex with volume: 

62(61 — Oi ) 2 (62 — a 2 )(b 3 — o 3 ) 2 /24. 

The polytope conv(T' U {w 7 }) = convj-u 1 , v 5 , v 6 , v 8 , v‘ } is a 4-simplex with volume: 

61(61 — tti) (62 — a 2 ) 2 ( 6 3 — a 3 ) 2 /24. 

To compute the volume of Vh, we sum the volume of the appropriate eight simplices, and we 
obtain the volume of Vh as stated in Theorem |4.1[ □ 
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Table 1 


Ineq 

Mivlvl 0 ) 

M{v 3 , W 3 1 ) 




17 


immediate 

immediate 

immediate 

immediate 

immediate 

is 


immediate 

immediate 

immediate 

immediate 

immediate 

i'l 


immediate 

by Lemma 

10.1 

immediate 

by Lemma 10.1 

by Lemma 

10.1 

21) 


immediate 

by Lemma 

10.1 

by Lemma 10.1 

immediate 

by Lemma 

11 i.l 

21 


immediate 

by Lemma 

lo.l 

immediate 

by Lemma 10.1 

by Lemma 

10.1 

22 


immediate 

by Lemma 

lo.l 

by Lemma 10.1 

immediate 

by Lemma 

li i.l 

2:1 


immediate 

immedia 

te 

immediate 

immediate 

immedia 

te 

2l 


immediate 

immediate 

immediate 

immediate 

immediate 

25 

26 

27 


see 

see 

imme 

35 

38 

dia 

te 

see i 10 . 

see 39 
immedia 

2 

te 

immediate 

see 36 
immediate 

see 36 
immediate 
immediate 

see [37| 

See Jla3 
immediate 

2s 


immediate 

immediate 

immediate 

immediate 

immediate 

23 


immediate 

immediate 

immediate 

immediate 

immediate 

:»l 


immediate 

immediate 

immediate 

immediate 

immediate 

31 


immediate 

immediate 

immediate 

immediate 

immediate 

2,2 


immediate 

immediate 

immediate 

immediate 

immediate 

2,2 


immediate 

immediate 

immediate 

immediate 

immediate 

2,1 


immediate 

immediate 

immediate 

immediate 

immediate 


6. Proof of Thm. 4A\ We compute the volume of the convex hull of the 12 extreme points 


which we claim are exactly the extreme points of system 3. In computing the volume of this 
polytope, we also prove that these are the correct extreme points and that we have therefore 
computed the volume of V 3 . 

The relevant points are the eight extreme points of Vh, plus an additional four points. Because 
we have already computed the volume of Vh , to compute the volume of V 3 , we need to compute 
the additional volume, compared with Vh, added by these four extra extreme points. To show that 
this is indeed the volume of V 3 , we keep track of which facets need to be deleted and added to the 
system of inequalities as we go. When this is complete, we have exactly system 3, and therefore we 
must also have the correct extreme points. 

We begin with system H from 12.2 As discussed in ij5j it would be convenient to add the four 
points to Vh separately; i.e., compute the additional volume each produces when added to Vh, and 
sum the results. To show that we can add two points separately and obtain the correct volume, we 
show that the intersection of the line segment between these points and Vh is non-empty (Lemma 


5.1). 


We show that we can add separately, ng 0 separately, and then v g 1 and vl 2 together by con¬ 
sidering the midpoints of the line segments between the relevant points. We consider L(vl,vl°), 
L(y fjUg 1 ), L(vfv}f). and L(y]f,v\ 2 ). We show that the midpoint of each line segment 

satisfies each of the inequalities of Vh by substituting this point into each inequality and checking 
the result. See Table [I] for a summary of the resulting substitutions. The table notes whether non¬ 
negativity of the resulting quantity follows immediately (after factoring), or by use of a technical 
lemma (after further explanation in the appendix), or after being rewritten in the way referenced 
in Figure [3} Because we have shown that each resulting quantity is nonnegative, we conclude that 
each of the midpoints intersect Vh, and therefore we can add uf separately, v\° separately, and 
then and Ug 2 together. 


6.1. Computing the (additional) volume of conv(V U {n®}). We now compute the 
additional volume of conv(V H U {vf}) compared to the volume of Vh- To do this we sum the 
volumes of conv({ng}UT) for each facet, F, of Vh such that n® is beyond that facet. We substitute 
vl into each inequality of system H , and we immediately see that it satisfies every inequality except 
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Figure 3 . For Table [l] 


(62 - a 2 )(bi - <u) (6163(i>2 - a 2 ) + a 2 a 3 (bi - ai)) 

2(6i 6 2 — aia 2 ) 

( 35 ) 

(b 2 b 3 - a 2 a 3 )(foi - ai) + (bib 3 - axa 3 )(fo 2 - a 2 ) 

2 

( 36 ) 

(62 -a 2 )(&i — Ri) (&i( 6 2 6 3 -a 2 o 3 ) + a 2 (6i& 3 -oia 3 )) 

2(6i b 2 — ai a 2 ) 

( 37 ) 

(b 2 - o 2 )(6i - a 1 )(ft 2 b 3 (6i - ai) + aia 3 (b 2 - a 2 )) 

2(6i b 2 — aio 2 ) 

( 38 ) 

(b 2 - a 2 )(bi - Oi) {b 2 (bib 3 - aio 3 ) + a\{b 2 b 3 - a 2 a 3 )) 

2(6i b 2 — ai a 2 ) 

( 39 ) 


From this we know that tif is beyond only one facet. The extreme points that lie on this facet 
are points v 1 ,v 2 ,v 6 and u 8 . The polytope conv{ti 1 ,ti 2 ,t’ 6 ,ti 8 ,tig} is a 4-simplex with volume: 

6102(61 — Oi) 2 (6 2 — a 2 ) 2 (63 — O3) 2 / (24(6i6 2 — Oia 2 )) . 

The facets of conv("P ff U{tif}) are the facets of Vh except inequality [25} We see this by computing 
the four additional facets that come from adding and noting they are already contained in 
system H: 


The facet through points v 1 , v 2 , v e and of is 

The facet through points v 1 , v 2 , v 8 and ti° is 

The facet through points v 1 , v 6 , v 8 and ti° is 

The facet through points v 2 , v 6 , v 8 and of is 


23 


31 


30 


27 


6.2. Computing the (additional) volume of conv(V U {u^}). We now compute the 
additional volume of conv('P// U {tig 0 }) compared to the volume of Vh- To do this we sum the 
volumes of conv({ti}°}U.F) for each facet, F, of Vh such that ti{° is beyond that facet. We substitute 
Ug° into each inequality of system H, and we immediately see that every inequality is satisfied 


except P?7> 


From this we know that ti{° is beyond only one facet. The extreme points that lie on this facet 
are points v 2 ,v 4 ,v 5 and v 6 . The polytope conv{ti 2 ,ti 4 ,ti 5 ,ti 6 ,ti{ 0 } is a 4-simplex with volume: 


ai6 2 (6i - ai) 2 (6 2 - a 2 ) 2 (6 3 - a 3 ) 2 /(24(6 i 6 2 - aia 2 )). 


The facets of conv('P ff U{ti} 0 }) are the facets of Vh except inequality [26} We see this by computing 
the four additional facets that come from adding tig 0 and noting that they are already contained 
in system H: 

• The facet through points v 2 , v 4 , v 5 and tig 0 is 

• The facet through points v 2 , v 4 , v 6 and vf is 

• The facet through points v 2 , v 5 , v 6 and tig 0 is 

• The facet through points u 4 , v 5 , v 6 and tig 0 is 

6.3. Computing the (additional) volume of conv('P^ U {tig } U {tig }). We now 

compute the additional volume of conv(T’ ff U{tig 1 }U{ti 3 2 }) compared to the volume of Vh- Because 
L(tig 1 ,tig 2 ) lies entirely outside of V3, we need to add them sequentially. 

We first compute the additional volume of con v(Vh U { ^3 1 }) compared to the volume of Vh- As 
we have done previously, we sum the volumes of convdti} 1 } U F ) for each facet, F, of Vh such 
that tig 1 is beyond that facet. We substitute tig 1 into each relevant inequality, and if the result is 


29 

24 

28 

32 
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negative then v\ 7 lies beyond that facet. It is immediate that v\ x violates inequalities 19-22 and 


satisfies inequalities [I7lfl8l [23p4| and [26p4| To see that inequality [25] is also satisfied see §10.4 


Therefore we have that v g 1 is beyond four facets, and we need to compute the volume of the 


convex hulls of U 3 1 with each of these facets. 


The extreme points that lie on the first facet are points u 1 ,?; 3 ,?; 5 and v 8 
convju 1 , v 3 , v 5 , v 8 , Ug 1 } is a 4-simplex with volume: 

aAOi - ai)(b 2 - a 2 ) 3 (b 3 - a 3 ) 2 /(24(M 2 - Oia 2 )) • 


The polytope 


The extreme points that lie on the second facet are points w 1 ,?; 4 ,?; 5 and v 7 . The polytope 
convty 1 ,v 4 ,v 5 ,v 7 ,vl 7 } is a 4-simplex with volume: 

ai& 2 (&i - oi ) 2 {b 2 - a 2 ) 2 (b 3 - a 3 ) 2 /{24^^ - a 3 a 2 )). 

The extreme points that lie on the third facet are points u 1 ,u 5 ,u 7 and v 8 . The polytope 
conv{u 1 ,u 5 ,'i; 7 ,'U 8 ,?; 3 1 } is a 4-simplex with volume: 

a ibi(b 3 — a 3 )(b 2 — a 2 ) 3 (b 3 — a 3 ) 2 / (24 (bib 2 — a 3 a 2 )). 

The extreme points that lie on the fourth facet are points u 1 ,^ 3 ,?; 4 and v 5 . The polytope 
convju 1 , v 3 , v 4 , v 5 , Ug 1 } is a 4-simplex with volume: 

ai5 2 (6i - a 3 ) 2 (b 2 - a 2 ) 2 (b 3 - a 3 ) 2 /(24^2 - aia 2 )). 

We now have a new polytope which is con v(V H U {ug 1 }). We refer to this polytope as T 3 , and we 
compute the facets of T 3 . 

We begin with the facets of Vh and delete the four facets that violated ( 19p2 ). Let us call 
this system T 3 ■ Now consider the four simplices we dealt with when computing the additional 
volume produced with Ug 1 . Each of these simplices has 5 facets; one of which corresponds to a 
deleted facet of Vh- 


The remaining 4 facets of the first simplex are described by the planes through the following 
sets of points: {u 1 ,^ 3 ,?; 5 ,^ 1 }, {u 1 , v 3 , v 8 , fg 1 }, {u 1 , v 5 , v 8 , Ug 1 } 
and {u 3 , v 5 , v 8 , fg 1 }. 

The remaining 4 facets of the second simplex are described by the planes through the following 
sets of points: {u 1 ,u 4 ,u 5 ,Ug 1 }, {u 1 , u 4 , v 7 , fg 1 }, {u 1 , v 5 , v 7 , Ug 1 } 
and {u 4 , v 5 , v 7 , fg 1 }. 

The remaining 4 facets of the third simplex are described by the planes through the following 
sets of points: {v 1 ,v 5 ,v 7 ,v\ 7 }, {u 1 , v 5 , v 8 , fg 1 }, {u 1 , v 7 , v 8 , Ug 1 } 
and {u 5 ,^ 7 ,?; 8 ,^ 1 }. 

The remaining 4 facets of the fourth simplex are described by the planes through the following 
sets of points: {v 1 : v 8 , {u 1 , v 3 , v 5 , fg 1 }, {u 1 , v 4 , v 5 , Ug 1 } 
and {u 3 ,^ 4 ,^ 5 ,^ 1 }. 

Consider these sixteen facets and exclude the facets that are shared by more than one simplex. 
This leaves eight facets. 

We can compute these eight facets to obtain the following: 

• The facet through points v 1 , v 3 , v 8 and Ug 1 is 

-— -(— a\ofb 3 + a\a 2 b 3 x 2 — a,ibi 0 ^ci 3 + ciib\ 0 %x 3 

b\b 2 — aia 2 

-\-CLiO, 2 b 3 Xi + ai6ia 2 6 2 a 3 + a,ibia 2 a 3 x 2 — aibici 2 b 3 x 2 — aibib 2 a 3 x 2 
T b 2 a 2 b 2 b 3 — b 2 a 2 b 2 x 3 — bia 2 b 2 b 3 x 1 — aiQ 2 f + b\b 2 f} > 0. 


(40) 
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The facet through points v 3 , v 5 , v 8 and v g 1 is 

/ - a 2 b 3 x 1 - aib 3 x 2 - bib 2 x 3 + a 3 a 2 b 3 + &i6 2 6 3 > 0. 

The facet through points v 1 , v 4 , v 7 and U 3 1 is 

/ - b 2 a 3 x 1 - b x a 3 x 2 - aia 2 x 3 + aia 2 a 3 + bib 2 a 3 > 0. 


(41) 


(42) 


,A />»H 

The facet through points v 1 , v 7 , v 8 and Ug 1 is 


The facet through points u , v 5 , v 7 and Ug 1 is 


32 


1 


(—Ui&iU 2 u 3 + (iibi(i\x 3 + ai&ia 2 a 3 x 2 — a,ibid 2 b 2 b 3 + a,iCi 2 b 2 b 3 Xi + bid 2 b 2 a 3 Xi 
+ b\a 2 b 2 b 3 — b\a 2 b 2 x 3 — b\Q, 2 b 2 b 3 Xi + b\b^a 3 — b\b 2 a 3 x 2 — bib^ci 3 Xi — CLid 2 f T b 3 b 2 f ) > 0. 


b 3 b 2 — a 1 a 2 


(43) 


The facet through points v 5 , v 7 , v 8 and Ug 1 is 
The facet through points v 1 , v 3 , v 4 and Ug 1 is 


The facet through points v 3 , v 4 


v 5 and Ug 1 is 


18 


17 


29 


There are four inequalities that are not already contained in system 7^ _ , we add these and in 
doing so obtain the system of inequalities that describes T 3 = conv(P ff U {ug 1 }). 

We now compute the additional volume of corw{T 3 U {ug 2 }) compared to the volume of T 3 - As we 


,12 


have done previously, we sum the volumes of conv({ug 2 } U F) for each facet, F, of F 3 such that v 
is beyond that facet. We substitute Ug 2 into each relevant inequality (i.e., the system of inequalities 
that describes %) and if the result is negative then Ug 2 lies beyond that facet. It is immediately 
clear that v\ 2 satisfies inequalities 


Ug 2 violates inequalities 40 and 43 


I7p8l [23p5l[27p4|and 


To see that inequality 


41 


26 


Therefore we see that v\ 2 is beyond two facets and we need to compute the volume of the convex 


42 We can also see immediately that 


is also satisfied see 110.4 


hull of Ug 2 with each of these facets. 


The extreme points that lie on the first facet are points v 4 ,v 3 ,v 8 and v\ 4 . The polytope 


convju 1 ,?; 3 ,'!; 8 ,^ 1 ,^ 2 } is a 4-simplex with volume: 


bia 2 (bi — ai) 2 (b 2 — a 2 ) 2 (b 3 — a 3 ) 2 / (24(bib 2 — a 3 a 2 )) 


The extreme points that lie on the second facet are points v l : v 7 ,v 8 and Ug 1 . The polytope 
conv{'u 1 ,u 7 ,u 8 ,U 3 1 ,Ug 2 } is a 4-simplex with volume: 

bia 2 (bi — ai) 2 ( 6 2 — a 2 ) 2 (b 3 — a 3 ) 2 / (24 {b\b 2 — aia 2 )). 

We now compute the additional facets; we take the four facets from adding each simplex and 
delete the facet that repeats. This leaves us with the following six facets to compute: 


1S 

v\ l and Ug 2 is 


42 


The facet through points v 1 , v 7 , v 8 and Ug 2 is 30 
The facet through points v 1 , v 7 , vl 1 and v\ 2 i 
The facet through points v 7 , v 8 
The facet through points v 1 , v 3 
The facet through points v 1 , v 3 


v 8 and v\ 2 is 
ul 1 and vl 2 is 


18 


31 


The facet through points v 3 , v 8 , Ug 1 and v\ 2 is 


17 


41 


By adding and deleting the appropriate facets to and from system H, we see that we arrive at 
system 3. 

Therefore, to compute the volume of V 3 , we sum the volume of Vh with that of the appropriate 
eight simplices, and we obtain our result. □ 
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Ineq 


M{v\ 


X°) 




Table 2 

M(vl,v\ 2 ) 


M{v\ 


X 1 ) 




20 


T? 




21 




24 


w> 


UU 


mi 


TTT 


7Z 


UU 


M 


immediate 

immediate 


see i 10.5 


immediate 
immediate 
see 


immediate 

immediate 


immediate 

immediate 


immediate 

immediate 


44 

and 

immediate 

see 

45 

and 

see ( 

10.11 


10.1 


see t 10.6 


Lemma 
see fiUT 


see 


Lemma 


45 


and 


Lemma 10.1 
immediate 


see 110.17 


see S 10.18 


see S10.8 


immediate 


see Slin 
immediate 


by Lemma 


46 


and 


10.1 



Lemma 

10.1 


TUT 


immediate 


by Lemma 


TUT 


47 


see 


immediate 

immediate 


48 


see 


27 


see 

50 


See ( 

10.10 


immediate 

see 

49 


see 

51 


immediate 

immediate 


49 


immediate 

immediate 

immediate 

immediate 

immediate 

immediate 


immediate 

immediate 


immediate 

immediate 

immediate 

immediate 

immediate 

immediate 


immediate 

immediate 

immediate 

immediate 


see i 10.13 


immediate 

immediate 

immediate 

immediate 

immediate 

immediate 


immediate 


see f 10.14 
immediate 
immediate 


immediate 

immediate 

immediate 

immediate 

immediate 

immediate 


immediate 

immediate 

immediate 

immediate 

immediate 

immediate 


7. Proof of Thm. 4.2 As with Theorem 4.4, we compute the volume of the convex hull of 


the 12 extreme points which we claim are exactly the extreme points of system 1. In computing 
the volume of this polytope, we also prove that these are the correct extreme points and that we 
have therefore computed the volume of V\. 

The relevant points are the eight extreme points of Vh, plus an additional four points. Because 
we have already computed the volume of Vh, to compute the volume of V± we need to compute 
the additional volume, compared with Vh, added by these four extra extreme points. To show that 
this is indeed the volume of Vi , we keep track of which facets need to be deleted and added to the 
system of inequalities as we go. When this is complete, we have exactly system 1 and therefore we 
must also have the correct extreme points. 


We begin with system H which can be found in 12.2 and we use the same principles as we used 
in the previous proof to compute the volume of V 3 . 

First we argue that we can add v\ to Vh separately, to Vh separately and then v\ l and 
v \ 2 together. To do this, we show that the midpoint of the line segment between v\ and all other 
additional points (u^°, v\ l and v\ 2 ) intersects Vh- We also show this is true for v\°. 

As in the previous proof we refer to the midpoint of the line between vl and v!? as M(vj , v *) 
and we show that the midpoint of each line satisfies each of the inequalities of Vh by substituting 
this point into each inequality and checking the result. See Table [2] for a summary of the resulting 
substitutions. The table notes whether nonnegativity of the resulting quantity follows immediately 
(after factoring), or by using a technical lemma, after further explanation in the appendix or after 
being rewritten in the way referenced in Figure [4j Because we have shown that each resulting 
quantity is nonnegative, we know that the midpoint intersects Vh, and therefore we can add u® 
separately, u(° separately, and then v\ x and together. 

7.1. Computing the (additional) volume of U }). We now compute the 

additional volume of con v(Vh U {u®}) compared to the volume of Vh- To do this, we sum the 
volumes of conv({u®}UT) for each facet, F, of Vh such that u® is beyond that facet. We substitute 


u® into the 18 relevant inequalities (17 34) and immediately see that it satisfies 17 19 23 26 and 


28 34 It is also immediate to see that inequality 27 is violated. To show that the remaining three 


inequalities are satisfied (20, 21 and 22) see 110.15 110.19 and 110.20 
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Figure 4 . For Table | 2 | 

(6 3 ~ a 3 )(6 2 - a 2 ) (6 3 (6ia 2 - ai6 2 ) + a 2 (&ia 3 - ai& 3 )) 
2(6263 — 0-20,3) 

(6i< 2 3 — ai6 3 )(6 2 — a 2 ) + (&ia 2 — ai6 2 )(6 3 — a 3 ) 

2 

(6 3 - a 3 )(6 2 - a 2 ) (6 2 (6ia 3 - gu 6 3 ) + a 3 (6ia 2 - ai& 2 )) 
2(6 2 6 3 — a 2 a 3 ) 

(6 3 - a 3 )(6 2 - a 2 ) (6j6 2 (6 3 - g 3 ) + aia 3 (6 2 - a 2 )) 
2(6 2 6 3 — a 2 a 3 ) 

(6 3 - a 3 )(6 2 - a 2 ) (6 2 (6 i 6 3 - aia 3 ) + a 3 (6i6 2 - aia 2 )) 
2(6 2 6 3 — a 2 a 3 ) 

(616 3 — cna 3 )(& 2 — a 2 ) + (616 2 — aia 2 )(& 3 — a 3 ) 

2 

( 6 3 - a 3 )( 6 2 - a 2 ) ( 6 j 6 3 ( 6 2 - a 2 ) + flia 2 ( 6 3 - g 3 )) 
2(6 2 6 3 — a 2 a 3 ) 

( 6 3 — a 3 )( 6 2 — a 2 ) ( 6 3 ( 6 i 6 2 - aia 2 ) + a 2 ( 6 i 6 3 -aia 3 )) 
2(6 2 6 3 — a 2 a 3 ) 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


From this we know that u® is beyond only one facet. The extreme points that lie on this facet 
are points v 2 ,v 3 ,v 6 and u 8 . The polytope conv{'f; 2 ,u 3 ,u 6 ,'U 8 ,u 5 ) } is a 4-simplex with volume: 

(^1 — Ol) 2 (^2 — 02) 2 (^3 — a 3 ) 2 /( 24 ( 62^3 — Gga 3 )) . 

The facets of conv((P ff Uju®}) are the facets of P# except inequality [27} We see this by computing 
the four additional facets that come from adding u® and noting that they are already contained in 
system H: 


The facet through points v 2 , 
The facet through points v 2 , 
The facet through points v 2 , 
The facet through points v 3 , 


v 6 and vf is 
v 8 and v\ is 
v 8 and vl is 
v 8 and vf is 


28 


31 


25 


34 


7.2. Computing the (additional) volume of conv(7U We now compute the 


additional volume of con v(Vh U {"c( 0 }) compared to the volume of Vr- To do this, we sum the 
volumes of conv({u(°}UF) for each facet, F, of Vh such that u 30 is beyond that facet. We substitute 


,10 


u t into the 18 relevant inequalities and immediately see that it satisfies |T7| , |T8| [20| [23] and 25p4 
is also immediate to see that inequality 


are satisfied (19 


It 

24 is violated. To show that the remaining three inequalities 


f 10.21 and 110.22 


21| and |22J) see §10.15 

From this we know that v\° is beyond only one facet. The extreme points that lie on this facet 


are points v 2 ,v 4 ,v 6 and v 7 . The polytope conv{u 2 ,ir 


,v 6 ,v 7 , 


Ug 0 } is a 4-simplex with volume: 


&2Q3(&1 — CLi) 2 (b -2 — G2) 2 (^3 — a 3 ) 2 /(24(62^3 — O- 2 CI 3 )) . 

The facets of couv(Vh U{u(°}) are the facets of Vr except inequality [24} We see this by computing 


v , v 


26 


the four additional facets that come from adding and noting that they are already contained 
in system H: 

• The facet through points v 2 , 

• The facet through points v 2 , 

• The facet through points v 2 , 

• The facet through points u 4 , v 8 , v 7 and v 4U is 


and v\° is 
v 7 and vl° is 
v 7 and vl° is 


33 


23 


32 
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7.3. Computing the (additional) volume of conv(V jj U U We now 

compute the additional volume of couv{Vh U{u( 1 }U{ , u( 2 }) compared to the volume of Vh- Because 
L(u( 1 ,u( 2 ) lies entirely outside of Vi, we need to add them sequentially. 

We first compute the additional volume of conv(V H U {u^ 1 }) compared to the volume of Vh- As 
we have done previously, we sum the volumes of convdu^ 1 } U F) for each facet, F, of Vh such 

into each relevant inequality and if the result is 
n 


,12 


that v \ 4 is beyond that facet. We substitute v 


negative then v \ 4 lies beyond that facet. It is immediate that v \ 4 satisfies inequalities [l7| [Tb| [23p6| 


and 28p4 In 110.16 we show that 27 is also satisfied. It is also immediate that v\ 4 violates the 
three facets described by 2p|22 We compute the volume of the convex hulls of with each of 
these facets. 

The extreme points that lie on the first facet are points v 4 ,v 4 ,v 3 and v 7 . The polytope 
conv{u 1 ,?; 4 ,?; 5 ,'U 7 ,'d 1 } is a 4-simplex with volume: 

b 2 d 3 (bi — di) 2 {b 2 — d 2 ) 2 (b 3 — d 3 ) 2 / (24(b 2 b 3 — d 2 d 3 )). 

The extreme points that lie on the second facet are points v 4 ,v 5 ,v 7 and v 8 . The polytope 
conv{w 1 ,u 5 ,T; 7 ,u 8 ,u( 1 } is a 4-simplex with volume: 

bia 3 (bx — ai)(i ) 2 — d 2 ) 3 (b 3 — d 3 ) 2 / {2A(b 2 b 3 — a 2 a 3 )). 

The extreme points that lie on the third facet are points v 1 ,v 3 ,v 4 and v 3 . The polytope 
conv{u 1 ,u 3 ,u 4 ,u 5 ,u( 1 } is a 4-simplex with volume: 

0162(61 “ 0 , 1 )(b 2 - a 2 ) 2 {b 3 - a 3 f /(24(6 2 6 3 - a 2 a 3 )). 

Unlike in system 3, we see immediately that there exists a fourth facet (described by |19[) which, 
under certain circumstances, v\ l is beyond. In particular, this is true if and only if cq b 2 b 3 — bia 2 a 3 > 
0. Therefore we continue with two cases. 


7.3.1. Case 1: ai&2^3 — ^l a 2 a 3 > 0 I n this case there exists a fourth facet (described by 
19) such that v \ 4 is beyond this facet. The extreme points that lie on this fourth facet are points 
v 1 ,v 3 , v 5 and v 8 . The polytope convju 1 , v 3 , v 5 , v 8 , uj 1 } is a 4-simplex with volume: 

{d\b 2 b 3 — b 1 a 2 a 3 )(bi — di)(b 2 — d 2 f{b 3 — a 3 ) 2 / (2A(b 2 b 3 — d 2 d 3 )). 

We now have a new polytope which is con v(Vh U (d 1 }) (in case 1). We refer to this polytope as 
Ti , and compute the facets of T 4 ■ 

We begin with the facets of Vh and delete the four facets that v \ 4 lies beyond (19 22). Let us call 
this system T 3 ~- Now consider the four simplices we dealt with when computing the additional 
volume produced with v\ 4 . Each of these simplices has 5 facets; one of which corresponds to a 
deleted facet of Vh- 


The remaining 4 facets of the first simplex are described by the planes through the following 
sets of points: {u 1 ,u 4 ,u 5 ,u( 1 }, {u 1 ,u 4 ,u 7 ,u( 1 }, {u 1 , v 5 , v 7 , v\ 4 } 
and {u 4 , v 5 , v 7 , uj 1 }. 

The remaining 4 facets of the second simplex are described by the planes through the following 
sets of points: {u 1 ,u 5 ,u 7 ,u( 1 }, {u 1 , v 5 , v 8 , uf 1 }, {u 1 , v 7 , v 8 , v\ 4 } 
and {v 3 ,v 7 ,v 8 ,v\ 4 }. 

The remaining 4 facets of the third simplex are described by the planes through the following 
sets of points: {u 1 ,u 3 ,u 4 ,u( 1 }, {u 1 , v 3 , v 5 , uj 1 }, {v 1 ^v 4 ^ 3 ,v\ 1 } 
and {u 3 , v 4 , v 3 , u} 1 }. 
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The remaining 4 facets of the fourth simplex are described by the planes through the following 
sets of points: {upupupup}, {upupuPup}, {upu 5 ,u 8 ,up} 
and {-u 3 , -u 5 , -u 8 , 


Consider these sixteen facets and exclude the facets that are shared by more than one simplex. 
This leaves eight facets. 

We compute these eight facets to obtain the following: 

• The facet through points u 1 , u 3 , v 8 and up is 


( — Oifl2®3^3 T O1G262O363 4 " & l < l 2 ( l 3 b 3 X 2 — O162O363X2 — + 610^03X3 


6 2 6 3 - a 2 a 3 v 1 1 1 — (52) 

+ 020363X1 + 610203X2 — b\a 2 a 3 b 3 x 2 -\-bia 2 b 2 b 3 — bia 2 b 2 b 3 x 3 — a 2 b 2 b 3 x 3 — a 2 a 3 f -\-b 2 b 3 f'j > 0 . 

• The facet through points u 3 , u 5 , v 8 and up is 


(—01020363 + a ,\( i 2 ( i 3 b 3 x 2 + ci \( i 2 b 2 b 3 x 3 + 01636 


b 2 b 3 — a 2 a 3 

— 016363X3 — O16263X2 + 0^0363X1 — b \ Ci 2 b 2 CL 3 b 3 + 61G262O3X3 
+ 63026263 — 61026263X3 — O26263X1 — CL 2 CL 3 f + 6263/) > 0. 


(53) 


The facet through points u 1 , v 4 , u 7 and up is 33 

The facet through points u 4 , u 5 , v 7 and up is 32 

The facet through points u 1 , v 7 , v 8 and v\ l is 

/ - 6263X1 - 6i0 3 x 2 - 6io 2 x 3 + b 3 a 2 a 3 + 6 i6 2 6 3 > 0. 

The facet through points v 5 , v 7 , v 8 and vl 1 is 


(54) 


The facet through points v 1 , v 3 , v 4 and vl 1 is 
The facet through points v 3 , v 4 


18 


17 


v 5 and v\ 4 is 


/ - a 2 a 3 xi - ai6 3 x 2 - ai6 2 x 3 + 010203 + ai6 2 6 3 > 0. 


(55) 


There are four inequalities that are not already contained in system T 3 ~, we add these and in 
doing so obtain the system of inequalities that describes T' = conv(T’ ff U {uj 1 }) (in case 1). 

We now compute the additional volume of conv(T ( 1 U {up}) compared to the volume of 7} 1 . As 
we have done previously, we sum the volumes of conv({u} 2 } U F) for each facet, F, of 7 } 1 such 
that up is beyond that facet. We substitute u ] 2 into each relevant inequality (i.e., the system that 
describes 7} 1 ) and if the result is negative then u ] 2 lies beyond that facet. It is immediately clear 
that u ] 2 satisfies inequalities [Tt| [l 8 | [23| [25p4| and 54p5 . We also see immediately that inequalities 
52 and 53 are violated. To see that inequality [24] is also satisfied see appendix §10.16 


Therefore we see that up is beyond two facets and we need to compute the volume of the convex 
hull of up with each of these facets. 

The extreme points that lie on the first facet are points v 4 ,v 3 ,v 8 and up. 


The polytope 


conv{u 1 ,u 3 ,u 8 ,up,u} 2 } is a 4-simplex with volume: 


0263(61 — Oi ) 2 (62 — a 2 ) 2 ( b 3 — a 3 ) 2 / (24(6263 — a 2 a 3 )) 


The extreme points that lie on the second facet are points u 3 ,u 5 ,u 8 and up. The polytope 
conv{u 3 ,u 5 ,u 8 ,up,u} 2 } is a 4-simplex with volume: 

0263(63 — Oi) 2 (6 2 — a 2 ) 2 ( b 3 — a 3 ) 2 /(24(6 2 6 3 — a 2 a 3 )). 


We now compute the additional facets; we take the four facets from adding each simplex and 
delete the facet that repeats. This leaves us with the following six facet defining inequalities to 
compute: 
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,11 


The facet through points v 1 , v 3 , v 8 and v \ 2 is 31 
The facet through points v 1 , v 
The facet through points v 1 , V 


and v } 2 is 


and v \ 2 is 


The facet through points v 3 , 
The facet through points v 3 , 
The facet through points v 5 , 


v 8 and v , 2 is 34 


v , v 


v\ l and v ] 2 is 


17 


54 


and v\" is 


55 


18 


By adding and deleting the appropriate facets from system H we see that we arrive at system 1. 
Therefore, to compute the volume of Vi, we sum the volume of Vh with that of the appropriate 
eight simplices, and we obtain our result for case 1 . 


7.3.2. Case 2: ai&2^3 — ^l a 2 a 3 — ^ I n this case it is immediate to see that v \ 4 satisfies 
[T9| and therefore lies beyond no further facets. This means we now have a new polytope which is 
con v(V H U {ug 1 }) (in case 2). We refer to this polytope as If, and we compute the facets of if. 


We begin with the facets of Vh and delete the three facets that v\ x lies beyond (20 22). Let us 


call this system lf~ . Now consider the four simplices we dealt with when computing the additional 
volume produced with vf. Each of these simplices has 5 facets; one of which corresponds to a 
deleted facet of Vh- 


The remaining 4 facets of the first simplex are described by the planes through the following 
sets of points: {v 1 ,v 4 ,v 5 ,vf}, {v 1 ,v 4 ,v 7 ,v\ 4 }, {v 1 , v 5 , v 7 , vf} 
and {v 4 ,v 3 ,v 7 ,vf}. 

The remaining 4 facets of the second simplex are described by the planes through the following 
sets of points: {u 1 ,v 5 ,v 7 {u 1 , v 5 , v 8 , uf 1 }, {v 1 ,v 7 ,v 8 ,v\ 4 } 

and {v 3 ,v 7 ,v 8 ,v\ 4 }. 

The remaining 4 facets of the third simplex are described by the planes through the following 
sets of points: {u 1 ,v 3 ,v 4 ,v \ 4 }, {u 1 , v 3 , v 5 , uf 1 }, {u 1 , v 4 , v 5 , uj 1 } 
and {u 3 , v 4 , v 5 , u} 1 }. 

Consider these twelve facets and exclude the facets that are shared by more than one simplex. 
This leaves eight facets. 

We compute these eight facets to obtain the following: 


• The facet through points v 1 , v 5 , v 8 and vf is 


• The facet through points v 1 , v 4 , v 7 and vf is 33 

• The facet through points v 4 , v 5 , v 7 and vf is 32 


-——-, (—01610503 T 01630203 X 2 + aibia 2 b 2 x 3 — 01636263 + 036363 X 1 + 630303 X 1 

62(61 ~ d \) 

— b 2 CL 2 62 O 3 — b 2 0 , 2 ( 13 X 2 — 61 O 262 O 3 X 1 T 6 j 0,2 6263 — b 2 O 262 X 3 — 63 O 2 62 63 X 3 — U 1 & 2 / + 63 62 /) ^ 0. 


(56) 


• The facet through points v 1 , v 7 , v 8 and vf 

• The facet through points v 5 , v 7 , v 8 and v \ 4 

• The facet through points v 1 , v 3 , v 4 and vf 

• The facet through points v 1 , v 3 , v 3 and vf 


is 

is 

is 

is 


54 

18 

17 


1 

a 3 (bi - oi 

+ 03630203 T 010 : 20363 X 1 — £ 11610363 X 2 + O 162 O 363 X 1 + 01616263 

- ££1616263X3 - 016263X1 - 610203X1 - Oia 3 / + 6 i0 3 /) > 0. 


) 


(—01020363 — O162O363 + 0^0363X2 + 0^6263X3 


(57) 
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The facet through points v 3 , v 4 , v 5 and vf is 55 


There are four inequalities that are not already contained in system T 2 ~\ we add these and in 
doing this we obtain the system of inequalities that describes T 2 = con v(Vh U {uj 1 }) (in case 2). 

We now compute the additional volume of con v(T[ 2 U {u( 2 }) compared to the volume of If. As 
we have done previously, we sum the volumes of conv({u( 2 } U F) for each facet, F, of Tf such 
that v\ 2 is beyond that facet. We substitute vf into each relevant inequality (i.e., the system of 
inequalities that describes T 2 in case 2) and if the result is negative then vf lies beyond that 
facet. It is immediately clear that vf satisfies inequalities 17, 18, 23, [25] 34, 54 and 55 We also see 


immediately that v\ 2 violates inequalities 19, 56 and 57 To see that 24 is also satisfied see { 10.16 


Therefore we know that v\ 2 is beyond three facets, and we need to compute the volume of the 
convex hull of vf with each of these facets. 

The extreme points that lie on the first facet are points u 1 ,?; 3 ,?; 5 and v 8 . The polytope 
conv{'u 1 ,u 3 ,u 5 ,u 8 ,u( 2 } is a 4-simplex with volume: 

(^1 — Q-l) 2 (^2 — 02) 2 (^3 — a 3 )7(24(M 3 — 0203 )) ■ 

The extreme points that lie on the second facet are points v 1 ^ 5 ^ 8 and v\ x . The polytope 
conv{'u 1 ,u 5 ,u 8 ,u( 1 ,w( 2 } is a 4-simplex with volume: 

(i>ia 2 (i>i — oi)(&2 — CL 2 ) 2 {b 3 — a 3 ) 3 /(24(& 2 6 3 — 0203 )) • 

The extreme points that lie on the third and final facet are points v 1 , v 3 , v 5 and vf. The polytope 
conv{'U 1 ,'U 3 ,u 5 ,'i; 41 ,' y i 2 } is a 4-simplex with volume: 

aib 3 (bi — Oi)( 6 2 — a2) 3 (^3 ~ a 3) 2 /(24(6 2 6 3 — o. 2 a 3 )) ■ 

We now compute the additional facets; we take the four facets from adding each simplex and 
delete the three facets that are repeated. This leaves us with the following six facet defining 
inequalities to compute: 


The facet through points v 1 , v 3 , v 8 and v} 2 i 
The facet through points v 3 , v 
The facet through points v 1 , v 


is 

v 8 and v\ 2 is 
11 


31 


34 


and v{ 2 is 

The facet through points v 5 , v 8 , v\ l and v{ 2 is 

The facet through points v 1 , v 3 , v\ 3 and v\ 2 is 

The facet through points v 3 , v 5 , vf and u( 2 is 


18 


17 


55 


tom system H we see that we also arrive at 


By adding and deleting the appropriate facets 
system 1 in case 2. 

Therefore, to compute the volume of V \, we sum the volume of Vh with that of the appropriate 
eight simplices, and we obtain our result for case 2. □ 


8. Proof of Thm. 4.3 
|4.3| immediately. □ 


A mapping from the proof of Theorem [T4] allows us to claim Theorem 


9. Possible extensions. Our results geometrically quantify the tradeoff between different 
convexifications of trilinear monomials. Of course it would be nice to use our results to develop 
guidelines for attacking trilinear monomials within an sBB code. In doing so, it should prove 
important to develop guidelines for how our results could be applied to formulations having many 
trilinear monomials overlapping on the same variables. We have seen that our results are very 
robust for scenarios where there is a high degree of overlap between trilinear monomials (see [27]). 
Also, we can imagine scoring each possible relaxation according to its volume, and then aggregating 
the scores to decide on what to do for each trilinear monomial. Another important issue is how 
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to effectively make branching decisions in the context of our relaxations. Guided by our volume 
results, we have made some significant progress in this direction (see mb¬ 
it would be natural and certainly difficult to extend our work to multilinear monomials having 
n > 3. In particular, advances for the important case of n = 4 could have immediate impact; [3] 
found, via experiments, that composing a trilinear and bilinear convexification in the manner 
suggested by ( x t Xj)xkXi was a good strategy. They further observed sensitivity to the bounds on the 
variables, but they reached no clear conclusion on how to factor in that aspect. Restricting to this 
type of convexification, we could apply our results by substituting w € [a,aj,b,bj\ to arrive at the 
trilinear monomial wXkXi, which can then be analyzed and relaxed according to our methodology. 
Of course, for a general quadrilinear monomial, there are six choices of which pair of variables will 
be treated as {xi,Xj}, so we can analyze all six possibilities and take the best overall. 

Also, there is the possibility of extending our results on trilinear monomials to (i) box domains 
that are not necessarily nonnegative, (ii) domains other than boxes, and (iii) other low-dimensional 
functions. 

We hope that our work is just a first step in using volume to better understand and mathemat¬ 
ically quantify the tradeoffs involved in developing sBB strategies for factorable formulations. 

10. Appendix. Throughout the proofs, we have repeatedly claimed that certain quantities 
are nonnegative for any choice of 01 , 02 , 03 , 61 , 62 , 63 , such that, 0 < Oj < b z , for all i and 


ai6 2 6 3 + 61O2O3 < 610263 + Oi6 2 a 3 < 6 i6 2 o 3 + 0i0 2 6 3 . 


In this appendix, we provide proofs for the cases that are not immediate. As will become apparent, 
we need to demonstrate that many different 6-variable polynomials are nonnegative on the relevant 
parameter space. Generally, such demonstrations can be tricky global-optimization problems, and 
in many cases sum-of-squares proofs are not available; rather, we often make somewhat ad hoc 
arguments. Still, we can place some efficiency on all of this by establishing some technical lemmas. 

10.1. We begin with the following lemmas that will be helpful in establishing the nonnega¬ 
tivity of certain quantities: 

Lemma 10.1. For all choices of parameters that meet our assumptions we have: b 1 a 2 — aib 2 > 0, 
b 1 a 3 — aib 3 > 0 and b 2 a 3 — a 2 b 3 > 0. 

Proof. (63 — a 3 )(b 3 ao — a 3 b 2 ) = bia 2 b 3 + a\b 2 a 3 — aib 2 b 3 — bia 2 a 3 > 0 by our original assumptions 
E This implies bia 2 — aib 2 > 0, because 63 — a 3 > 0. b 3 a 3 — a 3 b 3 > 0 and b 2 a 3 — a 2 b 3 > 0 follow from 
E in a similar way. □ 

Lemma 10.2. Let A , B, C,D,L,FgK with A > B > C > 0 and D > 0, E > 0, F < 0. Also let, 
D + E + F = 0. Then AD + BE + CF > 0. 

Proof. AD + BE + CF = AD + BE-C{D + E) > BD + BE - CD - CE = (B - C)(D + E) > 0. 

□ 

Lemma 10.3. Let A, B, C,D,E,F e M with A> B > C > 0 and D >0, E < 0, F < 0. Also let, 
D + E + F = 0. Then AD + BE + CF > 0. 

Proof. AD + BE + CF = —A(E + F) + BE + CF = E(B - A) + F{C - A) > 0. □ 

Lemma 10.4. Let A, B,C, D £ M with A> B > 0, C + D> 0, C > 0. Then AC + BD > 0. 

Proof. AC + BD>B{C + D)> 0. □ 
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10.2. Substituting M(v f,^ 1 ) into inequality 25 of the convex hull we obtain: 

(—0^02®3 "b ^ 1^26263 + Oifl2®3^1 "b 2aia2®3^1^2 — 2(2in2^1^2^3 

— O1O36162 — ®2®3^1 “1“ a 2^1^3 — fl2^1^2^3 + b^b^bfj ^2(6 i&2 — O1O2), 

the numerator of which can be rewritten as 

6163 ^(6163 — aia 3 )( b 2 — 02)^ + biO,2 ( fbiCL2 — O1&2X63 — °3)^ + O1O2 ^(6263 ~~ 0203)( a l — ; 


and is nonnegative by Lemmas 10.1 and 10.2 

10.3. Substituting M{v\ a . vlf) into inequality 26 of the convex hull we obtain: 

^ — Gqtt 2^3 ”1“ O2O2O362 — 0^0363 "b ^+6363 ~b 01036363 + 2O1O2O36162 
—20102616363 — 01616363 — 02036^62 + 6^6363) j c l { b \ b 2 — 010.2)5 

the numerator of which can be rewritten as 

6162 ^(6263 O2O3) (63 — Oi)) + 0362 ( fbiCl 2 — Oi& 2)(®3 — 63)^ + Oi02 ^(6363 — 0 i 03)(02 — 62)^ 5 


and is nonnegative by Lemmas 10.1 and 10.3 


10 . 4 . Substituting point u” into inequality 25 of the convex hull or substituting v\ l into 
inequality 26 we obtain: 

(— a\a\cLj, + 0^026263 + 01036363 + 80302036162 — 80302636263 
—O 1 O 361&2 — 02036^62 + 6 ^ 6363 ) j (6162 — O 1 O 2 ). 

The numerator can be rewritten as 

63 ^Oi 02 (oi &2 — 6162 ) + 6102(0302 — 0162 ) + 6162(6162 — 0303 )^ 

+ 03^0102(6162 — O1O2) + 6102(0162 — 6162) + 6162(0102 — 0162)^ =: 63T + 0,3Z. 

Then we can see X + Z = (6 2 — o 2 )(6i — Oi)(6i 6 2 — Oi0 2 ), which is nonnegative. Furthermore, by 
Lemma 10.3 we have X > 0. Therefore, by Lemma |10.4| the numerator is nonnegative. 


10.5. 


Substituting into inequality 19 of the convex hull we obtain: 

^ — 2 O 1 O 2 O 363 + O 1 O 263 — O 1 O 2 O 362 + 4 O 1 O 2 O 36263 — O 1 O 26263 — O 16263 + O 2 O 361 

+03036163 — 036363 — 40203616263 + 202616263 + 03616363) j 2(6363 — 0203), 

the numerator of which can be rewritten as 

6263 ^ 62(6303 — 0163 ) + 03(6163 — 26i0 3 + 0103 )^ + 0263 (fbi — Oi) (63 — a 2 )( 6 3 — 03 )^ 
+ O 2 O 3 ^ 62(20363 — 6163 — 0103 ) + 02(6103 — 0163 )^ =: 6263 X + 0263 !^ + (L 2 CI 3 Z. 


Now, we write X =: 6 2 X + a 2 lX and see that lb + IF = (61 — Oi)(6 3 — a 3 ) > 0 and, by Lemma 10.1 


V > 0. Therefore X > 0 by Lemma 10.4. Because X + X + Z = 0, by Lemma 10.2 we have that the 
numerator is nonnegative. 
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10.6. Substituting M(v\,v\°) into inequality |20| of the convex hull, we obtain: 

(—01020363 — 203020362 + 40302036263 + 030362 — 03036363 — 036363 + 030363 

+02036362 — 40203636263 + 02636363 — 036363 + 203636363) j 2(6363 — 0203), 

the numerator of which can be rewritten as 

6263 ^63(6302 — 0363) + 03(6362 — 26302 + 0302)^ + 62O3 ^(63 — O3) (63 — 02)(6 3 — 03)^ 

+ O 2 O 3 ^ 63(20363 — 6362 — 0302 ) + 03(6302 — 0362 )^ =• b 2 b 3 X + b 2 a 3 Y + a, 2 a 3 Z. 


Now, we write X =: 63 Y + a 3 lY and see that 14 + 1+ = (63 — Oi)( 6 2 — a 2 ) > 0 and, by Lemma 10.1 


V > 0. Therefore X > 0 by Lemma [10.4[ Because X + Y + Z = 0, by Lemma [10. 2 1 we have that the 
numerator is nonnegative. 


10.7. Substituting M(vl,v\ x ) into inequality 20 of the convex hull, we obtain: 


(—03030363 + 20302036263 — 030363 + 03036363 — 036363 + 030363 — 03036363 
—20203636363 + 02636263 + 036363) 72(6263 — 0303), 


the numerator of which can be rewritten as 


6263 (( 6 ia 2 - ai 6 2 )( 6 3 - a 3 )) + 6 2 a 3 (( 6 2 a 3 - a 2 b 3 )(bi - 03 )) + a 2 a 3 ((6303 - ai 6 3 )(a 2 - & 2 )) , 


which is nonnegative by Lemmas 10.2 and 10.1 


10.8. Substituting M^vfvl 1 ) into inequality 21 of the convex hull we obtain: 

(62 — 02) (03030363 — 036263 — 030363 — 02036363 + 026363 + 036363) j 2(6263 — 0303), 
where the second multiplicand of the numerator can be rewritten as 

63^63(63a 2 - 0362)) + 03(030263 - 630203 + 6i6 2 a 3 - 630263^) =: 6 3 Y + a 3 Z, 


now Y + Z = (63 a 3 — 0363 X 62 — a 2 ) > 0 (Lemma 10.1), and Y > 0 (Lemma 10.1), therefore by 
Lemma 10. 4| we have that 63 Y + a 3 Z is nonnegative. 


10.9. Substituting M(vfv\ l ) into inequality 22 of the convex hull we obtain: 


(63 — 03 ) (030363 — 03036363 + 030363 — 036363 — 0^0363 + 02636263 ) j 2(6363 — a 2 a 3 ), 
where the second multiplicand of the numerator can be rewritten as 

62 (^ 63(6302 — 0362 ) + 03(6303 — a 2 b 3 ) \ + 02 ^ 02(0363 — bia 3 )\ =: 6 2 Y + 02 ^, 


where Y + Z = (6302 — 0362 X 63 — a 3 ) > 0 and Y > 0 (both Le mm a 10.1), therefore by Lemma 10.4 
we have that this term is nonnegative. 
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10.10. Substituting into inequality [27] of the convex hull we obtain: 

(— O 1 O 2 O 3 4“ O 1 O 2 O 363 — 030363 4“ 2 O 3 O 2 O 36263 — O 3 O 36263 T ^ 2^1 ^3 


+020361^2 — 2 (X 2 (l 3 bib 2 b 3 — CL 2 bib 2 b 3 4 “ ^1X2^3) j 2(6363 — O2O3), 
the numerator of which can be rewritten as 

6263 ^(6163 — aia 3 )(6 2 — <*2)^ + o , 2 b 3 [{ b 2 a 3 — CL 2 b 3 )(ai ~ 61)^ + <^2^3 ^(6362 — aia 2 )(a 3 — 63)^ . 
which is nonnegative by Lemmas |1 0.3 and 10.1 


10.11. Substituting M(uJ°,uJ 2 ) into inequality 19 of the convex hull we obtain: 

(— O3O263 — O3O2O362 4 “ 2O3O2O36263 4 “ CL \ Oj 2 b 2 b 3 — O36963 -\- 0 , 2 Cl 3 b \ 

—03036363 + 036163 — 2 a 2 a 3 bib 2 b 3 + 03616363) j 2(6263 — 0203), 

the numerator of which can be rewritten as 

6263 ((61 a 3 - ai6 3 )(6 2 - a 2 )) 4 - a 2 6 3 ((6 2 a 3 - a 2 6 3 )(ai - 61)) + a 2 a 3 ((6ia 2 - ai6 2 )(a 3 - 63)). 
which is nonnegative by Lemma 10. 3| and Lemma |10.1| 

10.12. Substituting M{v\°,v\ 2 ) into inequality [ 2 T] of the convex hull we obtain: 

(63 — 03) (03020362 — 016363 — 030361 + 036363 — O2O36162 + 036163) j 2 ( b 2 b 3 — 0203), 

where the second multiplicand of the numerator can be rewritten as 

62 ^62(6303 — Oi 6 3 )^ + 02^63(0263 — 6203) + 03(0362 — 61 a 2 )^ =: 6 2 Y + o 2 Z, 


where Y + Z = ( 6 ia 2 — ai 6 2 )( 6 3 — o 3 ) > 0 and Y > 0 (both by Lemma 10.1). Therefore by Lemma 


10.4 we have that 6 2 Y 4 -a 2 Z is nonnegative. 

10.13. Substituting M(oJ°,uJ 2 ) into inequality 22 of the convex hull we obtain: 

(6 2 — o 2 ) (oi 0 2 6 2 + 030362 — 03036263 — 036263 — a 2 a 2 6i + 03636263) j 2(6263 — < i 2 cl 3 ), 

where the second multiplicand of the numerator can be rewritten as 

63 (cLi(a 2 b 3 — 6203 ) + 6 2 ( 6 i 0 3 — 0163 )^ + 03 ^ 03(0362 — bia 2 )) =: 63 Y + a 3 Z, 

where Y 4 - Z = (6303 — Oi6 3 )(6 2 — a 2 ) > 0 and Z < 0 =>• Y > 0 (both by Lemma 10.1). Therefore by 
Lemma 10.4 we have that 6 3 Y + a 3 Z is nonnegative. 

10.14. Substituting M(v\°,v\ 2 ) into inequality 24 of the convex hull we obtain: 

(—O1O2O3 T O1O2O362 T 201O2O36263 — CL \ Ci 2 b 2 b 3 — O1O362 4 ~ U2O36163 
—2 a 2 a 3 bib 2 b 3 4- 036163 — 03636363 4- 636363) j2 ( b 2 b 3 — a 2 a 3 ), 

the numerator of which simplifies to 

6263 ((6 i6 2 - 0302X63 - o 3 )) 4- 6 2 a 3 ^(6 2 a 3 - a 2 6 3 )(6i - 03)) + o 2 a 3 ((6363 - Oia 3 )(a2 - & 2 )) , 
and is nonnegative by Lemmas 10.1 and 10. 2| 
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10 . 15 . 


into inequality [l 9 | we obtain: 


Substituting point 7+ into inequality 20 of the convex hull or substituting point v\ 


10 


(—01020363 — 03020362 + 30102036263 — 016363 + a 2 2 a\bi 
—80203616263 + 02616263 + 03616363) j {b 2 b 3 — 0203), 


the numerator of which can be rewritten as 

63 ^62(62(6303 — 0163) + 02(0303 + 6363 — 26303))^ 

O3 ^2(62(20163 — 0i03 — 6363) + 02(6303 — 


+ ■ 


—: 63 Y + a^Z. 


10.1 


Now, we write Y =: 6 2 Y + a 2 b 2 W and see that V + W = (6 X — Oi)( 6 3 — o 3 ) > 0 and, by Lemma 
V > 0 . Therefore Y > 0 by Lemma 10.4 Because Y + Z = (61 a 3 — Oi 6 3 )(6 2 — o 2 ) 2 > 0 (Lemma 10 TJ), 
by Lemma 10.4 we have that the numerator is nonnegative. 


10 . 16 . 


Substituting point v 3 into inequality 27 of the convex hull or substituting point v\ 


12 


into [ 24 ] we obtain: 

(-oia^o 2 + 30103036263 - 03026263 — 01036363 + 03036163 
+02036363 — 30203636263 + 636363) ^(6263 — a 2 o 3 ), 

the numerator of which can be rewritten as 

61 ^6263(6263 — 0203) + 6203(0203 — 0263) + 0203(0263 — b 2 b 3 fj 

+ <+ ^6263(0203 — a 2 b 3 ) + 6303(0263 — 6263) + 0203(6263 — 0203)^ =: 63+ + Oi-Z, 


where: Y + Z = (6 2 6 3 — a 2 a 3 )(b 2 — a 2 ){b 3 — a 3 ) > 0 , and by Lemma 10.3 we have Y > 0 . Therefore 
by Lemma 10 . 4 | we have that the numerator is nonnegative. 


10 . 17 . Substituting M(y\,v\°) into inequality 21 of the convex hull we obtain 

(26263 — 0263 — 62O3) {^d\d 2 (x 3 — 016263 — 2 b\a 2 a 3 + b 3 a 2 b 3 + b\b 2 a 3 ) j 2(6263 — a 2 a 3 ), 
where the second multiplicand of the numerator can be rewritten as 


62(6103 — ciib 3 ) + 02(0303 — 26 i 0 3 + 6163), 


which is nonnegative by Lemmas 10.1 and 10.4 


10 . 18 . Substituting M(ui,u(°) into inequality 22 of the convex hull we obtain 

(■ a 2 b 3 + 6 2 a 3 - 2 a 2 a 3 ){aia 2 b 3 + ai6 2 a 3 - 2ai6 2 6 3 - b 3 a 2 a 3 + b 3 b 2 b 3 ) j 2(6 2 6 3 - o 2 o 3 ), 

where the second multiplicand of the numerator can be rewritten as 

62(6363 — 20163 + O1O3) + 02(0363 — 6303), 


which is nonnegative by Lemmas 10.1 and 10.4 
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10.19. 


Substituting point v\ into inequality 21 of the convex hull we obtain 


63(62 — a2) ^62(630.3 ~ 0-363) + 02(01123 — 26 i 0 3 + 6363)^ j (6263 — O2O3)) 

which is nonnegative by Lemmas 10.1| and |10.4[ 

10.20. Substituting point u® into inequality [22] of the convex hull we obtain 

^2(63 — O3) ^62(6363 — 20363 + OiOg) + 02(0363 — 6303)^ ^(6263 — 0203), 


which is nonnegative by Lemmas 10.1 and 10.4 


10.21. Substituting point into inequality 21 of the convex hull we obtain 

62(63 - 03)^62(6303 - 0363) + 02(0303 - 26303 + 6363)^ j (6263 - 0203), 


which is nonnegative by Lemmas 10.1 and 10.4 


10.22. Substituting point v{° into inequality 22 of the convex hull we obtain 

03 (62 — 02) ^62(6363 — 20363 + 0303) + 02(0363 — 6303)^ j (6263 — 0203), 


which is nonnegative by Lemmas 10.1 and 10.4 
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